Band Offsets at the Si/SiC>2 Interface from Many-Body Perturbation Theory 



OO 
O 

o 

(N 



(N 



R. Shaltaf, 1 G.-M. Rignanese, 1 X. Gonze, 1 Feliciano Giustino, 2, 3, B and Alfredo Pasquarello 2 ' 3 

'European Theoretical Spectroscopy Facility (ETSF) and Unite Physico-Chimie et de Physique des Materiaux (PCPM), 
Universite catholique de Louvain, 1 Place Croix du Sud, B-1348 Louvain-la-Neuve Belgique 
Ecole Polytechnique Federale de Lausanne (EPFL), 
Institute of Theoretical Physics, CH-1015 Lausanne, Switzerland 
Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), CH-1015 Lausanne, Switzerland 

(Dated: May 2, 2008) 

We use many-body perturbation theory, the state-of-the-art method for band gap calculations, to 
compute the band offsets at the Si/SiC>2 interface. We examine the adequacy of the usual approx- 
imations in this context. We show that (i) the separate treatment of band-structure and potential 
lineup contributions, the latter being evaluated within density-functional theory, is justified, (ii) 
most plasmon-pole models lead to inaccuracies in the absolute quasiparticle corrections, (iii) vertex 
corrections can be neglected, (iv) eigenenergy self-consistency is adequate. Our theoretical offsets 
agree with the experimental ones within 0.3 eV. 

PACS numbers: 71.10.-w,73.40.Ty,73.40.Lq,85.40.-e,85.60.-q 
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The band offsets (BOs) are among the most important 
properties of a heterostructure. Their precise knowledge 
is extremely important to engineer electronic [lj and op- 
toelectronic devices. Various theoretical methods have 
been used to estimate the BOs 0, 0, H|. Among these, 
density functional theory (DFT) calculations allow one to 
study the variation of interface dipoles with the interface 
structure, without resorting to experimental input 

In the DFT approach, the valence band offset (VBO) 
and the conduction band offset (CBO) are conveniently 
split into two terms: (V,C)BO = AE®J T + AV. The 
first term AE® FT > referred to as the band- structure con- 
tribution, is defined as the difference between the valence 
band maxima (VBM) or the conduction band minima 
(CBM) relative to the average of the electrostatic poten- 
tial in each material. These are obtained from two inde- 
pendent standard bulk calculations on the two materials. 
The second term AV, called the lineup of the average of 
the electrostatic potential across the interface, accounts 
for all the intrinsic interface effects. It is determined 
from a supercell calculation with a model interface. 

However, the DFT-BOs suffer from two important lim- 
itations, namely, the well-known DFT band-gap prob- 
lem, and the use of an approximate functional to model 
the exchange-correlation (XC) energy, such as the local 
density approximation (LDA) or the generalized gradi- 
ent approximation (GGA). These affect the value of the 
VBM and the CBM, and consequently the calculated 
BOs. Orbital-dependent approaches [|| and many-body 
perturbation-theory (MBPT) within the GW approxima- 
tion have been used to try to overcome these short- 
comings d, 0, [l(| ■ While the former may be less compu- 
tationally demanding, their reliability cannot be assessed 
a priori ;f| in contrast with MBPT. In Refs. [l,[Ioj|, the 
lineup of the potential AV was approximated by its DFT 
value, arguing that it depends primarily on the charge 
distribution at the interface which is a ground state quan- 



tity, hence little affected by many-body effects. Doing so, 
only the band-structure contribution is modified: 



A^f = A<r 



A(5E V , C 



(1) 



where 8E V (6E C ) is the quasiparticle (QP) correction at 
the VBM (CBM): 5E, = Ef p - Ef FT for i = v,c. It 
is important to stress that these corrections, which are 
obtained from bulk calculations, are the only additional 
ingredients that are required when DFT calculations of 
the VBO and CBO already exist. 

Interestingly, for various semiconductor interfaces, the 
QP corrections of the band edges are found to be almost 
the same on both sides [l(| leading to A(SE V ) < 0.2 eV 
in Eq. ([I]). As a result of this cancellation of errors, DFT 
is quite successful for the same interfaces [lH with errors 
ranging from 0.1 to 0.5 eV, despite its limitations men- 
tioned above. This relative success of DFT explains why 
it has been widely used to predict the VBO for a wide 
range of interfaces. And, when needed, the CBO was 
also predicted using a simple scissor operator to correct 
the band gap to the experimental value. This assump- 
tion was further motivated by the fact that MBPT calcu- 
lations going beyond GW by including an approximate 
vertex correction (GWT) showed that the valence-band- 
edge energy remained at its DFT value for silicon, the 
whole correction going to the conduction bands 0,111. 

However, when it comes to semiconductor-insulator or 
insulator-insulator interfaces, it appears that the DFT 
errors on the VBO can be much more important. For in- 
stance, for the Si/Si02 interface, the VBO are calculated 
to be 2.3-3.3 eV [!, 0, [H, [l6| in noticeable disagree- 
ment with the experimental results of 4.3 eV [13] • It 
seems that, for such systems, the cancellation of errors is 
not as good, emphasizing the need to go beyond DFT by 
including QP corrections. The same is probably true for 
interfaces found in transistors or organic photovoltaics. 
Before applying such a highly-demanding method in a 
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predictive way, its quality needs to be assessed for an 
interface for which the VBO and CBO are well-known 
experimentally, such as the Si/Si02 one. 

In this Letter, the effect of QP corrections on the VBO 
and CBO at the Si/Si02 interface is investigated. By per- 
forming a quasiparticle self-consistent GW (QSGW) cal- 
culation on a small interface model, we demonstrate that 
the line up potential is well described in DFT. For two 
more realistic models, the BOs calculated within DFT 
are corrected by computing the QP corrections for bulk 
Si and Si02. The latter are found to be crucial to reach 
a good agreement with experimental results. Their effect 
can not be accounted for by a simple scissor operator on 
top of DFT results. The inclusion of the vertex correc- 
tion is shown not to affect the GW results significantly, 
whereas QSGW does not improve the results. 

Three different interface models are adopted in which 
ideal /3-cristobalite is matched to the Si(100) surface, tak- 
ing the theoretical value of bulk Si (asi—5AS A) as the 
lattice constant. The models I and II were generated in 
a previous work by two of us in order to calculate 
the BOs within DFT. They consist of 11 Si monolay- 
ers and up to 10 Si02 molecular layers which guarantees 
well converged BOs. Model III is limited to 8 Si mono- 
layers and 8 Si02 molecular layers in order to perform a 
QSGW calculation which is quite demanding computa- 
tionally. For bulk Si, the theoretical value asi is adopted 
for the lattice constant. For Si02, we consider on the one 
hand the cubic structure (space group Fd3m) whose re- 
laxed lattice constant is asi0 2 =7.43 A; and, on the other 
hand, a tetragonal structure (space group IAi/amd) in 
which two sides are strained to match the Si lattice con- 
stant asio 2 —^ a Si as in the interface models of Ref. [l6j |. 
while the third is left free to relax reaching a value of 
c Si0 2 =6.59 A. In both structures, the Si-O-Si bond an- 
gle is 180°. Hereafter, the cubic and strained structures 
will be referred to as c-Si02 and s-Si02, respectively. 

All our calculations are performed using the ABINIT 
package jl^|. Only valence electrons are explicitly con- 
sidered using norm-conserving pseudopotentials [l9( to 
account for core-valence interactions. By coherence with 
the DFT calculations of the B Os fig , the XC energy is 
also described within the GGA [20[ . 

At the DFT level, the band gaps E® FT for the bulk 
systems are found to be 0.7 eV for Si, 5.4 eV for c-Si02, 
and 5.1 eV for s-Si02 [I]}. For the interface models I and 
II, the VBOs (CBOs) are calculated to be 2.6 and 2.5 eV 
(1.6 and 1.8 eV), respectively [22| . 

Three approaches are used to compute the QP cor- 
rections. In the first approach (GW), the self-energy 
is calculated self-consistently by updating the eigenen- 
ergies in both the dielectric matrices and the Green's 
functions while keeping the DFT wavefunctions. In the 
second approach (QSGW), both the eigenenergies and 
the wavefunctions arc updated [23j |. In both approaches, 
the dielectric matrices are computed in the random-phase 



approximation (RPA) using the sum over states formu 
lation 
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In the third approach, a vertex correction 
r is included in both the screening and the self-en ergy 
by going beyond the RPA and including XC effects [1 21 ] . 
and only the eigenergies are updated. Details about the 
DFT-GGA kernel are provided in Ref. 

Within GW, the frequency dependence of the dy- 
namically screened Coulomb potential W is most of- 
ten approximated u sing various plasmon pole models 

(ppms) r~ 



27, 2i 



|29| . The advantage is not only to 
reduce the computational load, but also to obtain an an- 
alytic expression for the self-energy. In fact, PPMs have 
proven to be very effective in producing band gaps in 
good agreement with experiments. 

However, while the QP correction to the gap 8E g 
was found not to be very sensitive to the choice of the 
PPM [3(| , we observe that SE V and 6E C may vary from 
one PPM to another, as reported in Table U for Si and c- 
Si02- In particular, the variation of 5E g with the PPM 
is more pronounced in c-Si02 (up to 0.4 eV difference 
between the extreme cases) than in Si. Since a precise 
knowledge of the QP corrections at the band edges is re- 
quired for the band offsets calculations, it is necessary to 
go beyond PPMs taking explicitly into account the fre- 
quency dependence of W. 



TABLE I: Quasiparticle corrections (in eV) at the VBM 
(6E V ), at the CBM (5E C ), and for the band gap (SE g ) for 
Si and c-Si02- The corrections are calculated within GW us- 
ing the PPMs proposed by Hybertsen and Louie (HL) [26j |. 
von der Linden and Horsch (vdLH) [27l |. Godby and Needs 
(GN) 0, Engel and Farid (EF) [3], and without PPM. 





HL 


vdLH 


GN 


EF 


no PPM 


Si SE V 


-0.6 


-0.6 


-0.4 


-0.6 


-0.4 


5E C 


+0.1 


+0.1 


+0.2 


+0.1 


+0.2 


SEg 


+0.7 


+0.7 


+0.6 


+0.7 


+0.6 


c-Si0 2 5E V 


-2.6 


-2.5 


-2.0 


-2.3 


-1.9 


SE C 


+1.3 


+1.1 


+1.5 


+ 1.2 


+1.5 


5E g 


+3.9 


+3.6 


+3.5 


+3.5 


+3.4 



In this work, the explicit frequency dependence is 
obtained using the deformed contour integration tech- 
nique (HI. The calculated QP corrections are reported 
in Table [TT] for Si, c-Si0 2 , and s-Si0 2 . 

The QP gaps can easily be obtained from these results 

E QP = E DFT + 5Eg _ 

and 8.8 eV for c-Si0 2 



Within GW, E$ p =1.3 eV for Si 



These results are in good agree- 
ment with the experimental values of 1.2 and 8.9 eV [32} , 
respectively. While GWT results (1.3 and 8.5 eV) are 
very close, the QSGW results (1.5 and 9.5 eV) are over- 
estimated like systematically observed in Ref. [33| ■ 

Turning to SE V and 5E C , the comparison of our results 
with previous calculations is not always straightforward. 
In most of the cases, the focus is on the band gap and 
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TABLE II: Quasiparticle corrections (in eV) at the VBM 
{6E V ), at the CBM {5E V ), and for the band gap (8E g ) for 
Si, c-Si02, and s-SiC>2. The corrections are calculated within 
GW, GWT, and QSGW. 



Si 



c-Si0 2 



s-Si0 2 



SE V 
SE C 
SE n 



GW GWT QSGW GW GWT QSGW GW GWT QSGW 



-0.4 +0.1 -0.6 

-0.2 +0.7 +0.2 
-0.6 +0.6 +0.8 



-1.9 -1.3 -2.8 
+1.5 +1.8 +1.3 
+3.4 +3.1 +4.1 



-1.9 
+1.4 
+3.3 



-1.3 -2.8 
-1.8 +1.1 
-3.1 +3.9 



the QP corrections are not given explicitly (the VBM is 
set to zero after correction). Besides, 8E V and 6E C are 
very sensitive to the degree of convergence reached (the 
QP gap converges much faster than the QP corrections) , 
not to mention the effect of the PPM (see Table Q} . For 
Si, our GW and GWT results compare quite well with 
those of Fleszar and Hanke [l3| which also do not rely 



on any PPM [34(. No such comparison can be made for 
c-SiC>2. Note however that the QP corrections for c-SiC>2 
and s-SiC>2 differ by up to 0.2 eV, indicating a dependence 
of SE V and SE C on the strain. This is at variance with 
the findings of Ref. [35[ in the case of Si under isotropic 
strain. A possible explanation for this difference may be 
that the strain is not isotropic in s-SiC>2 
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Finally, 

note that the PPM proposed by Godby and Needs [28[ 
leads to QP corrections (see Table |T]) in excellent agree- 
ment with those of the contour deformation method, at 
variance with the other PPMs. This finding might be 
generalized after proper investigation. 

While GW and QSGW lead to a lowering of the VBM 
of Si (slightly larger for QSGW) compared to the DFT 
result, the inclusion of the vertex correction brings it back 
to roughly its original value with a small shift of 0.1 eV, 
all the QP correction being on the conduction band. A 
similar result was also found previously [Til . [i~3T ] giving 
some motivation to the use of a scissor operator to com- 
pute the CBO within DFT. For SiC>2, our results are very 
different. First, the VBM is also raised when going from 
GW to GWT, but it definitely does not reach the DFT 
level back. This indicates the recovery of the DFT VBM 
with GWT is a coincidence in Si. And, it definitely rules 
out the use of a simple scissor operator for the compu- 
tation of the BOs. Second, the lowering of the VBM in 
QSGW in much larger than in GW (by about 0.9 eV). 
In fact, it seems that the systematical overestimation of 
the band gap in QSGW originates essentially from a too 
strong lowering of the VBM. Note that the increase of 
the CBM is slightly lower in QSGW than in GW. 

In order to compute the QP corrections to the BOs, 
the many-body effects on Ay also need to be investi- 
gated. This is done by comparing the electronic density 
and the resulting AV calculated within DFT and QSGW 
for model III. In Fig.[lja), the difference between the pla- 
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FIG. 



1: Difference between DFT and QSGW calculations for 
model III of (a) the planar average of the electronic density 
and (b) the macroscopic average of the local potential. The 
density is expressed in me/a.u., and the potential in meV. 



nar average of the DFT electronic density pdft(z) and 
the QSGW one pcw(z) is reported. The QSGW result 
differs only slightly from DFT, with a maximum change 
of about 1 me/a.u. in the interface region. The differ- 
ence between the QSGW and DFT macroscopic average 
of the local potential Vqw — Vdft ; which is reported in 
Fig. QIb), is about ~ 45 meV in the interface region, and 
less than 12 meV in the bulk regions. This gives rise to 
a net difference of ~ 20 meV in the lineup of the poten- 
tial AV. Such a small difference in AV suggests that 
the interfacial charge density and, consequently, the as- 
sociated dipole moments are well described within DFT. 
This justifies the assumption that the line up potential 
can be taken to be the same as in DFT. 

Finally, using Eq. |T]), the band offsets can be com- 
puted within MBPT at the GW, GWT, and QSGW lev- 
els. Our results arc reported in Table IIIII and compared 
with the experimental ones. Within GW the agreement 
is excellent for both the VBO and CBO (less than 0.3 
eV difference) . The effect of the vertex correction is less 
than 0.1 eV on the BOs. This results from a cancellation 
of the effects on each side of the interface. At this stage, 
we cannot say whether this result can be generalized to 
any interface. In contrast, the effect of quasiparticle self- 
consistency is more pronounced. Due to the large low- 
ering of the VBM of Si02, the VBO is increased by 0.7 
eV compared to GW, leading to an overestimation by up 
to 0.5 eV. At the same time, the CBO is slightly smaller 
than in GW leading to an underestimation by up to 0.6 
eV. So that QSGW is the worst of the three approaches. 

In summary, we have investigated the band offsets at 
the Si/Si02 interface using many-body perturbation the- 
ory. Starting from the BOs obtained within DFT for two 
model interfaces, the QP corrections have been computed 
within eigenenergy self-consistent GW and GWI 1 , and 
quasiparticle self-consistent GW taking the frequency de- 
pendence of the screened potential explicitly into account 
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TABLE III: Quasiparticle band offsets (eV) for cubic and 
strained SiQ 2 using GW, GWV, and QSGW. 











cubic 






strained 






Model 


DFT 


GW 


awr QSGW 


GW 


awr 


QSGW 


Expt. 


VBO 


I 


2.6 


4.1 


4.0 


4.8 


4.1 


4.0 


4.8 


4.3 




II 


2.5 


4.0 


3.9 


4.7 


4.0 


3.9 


4.7 




CBO 


I 


1.6 


2.9 


2.7 


2.7 


2.8 


2.7 


2.5 


3.1 




II 


1.8 


3.1 


2.9 


2.9 


3.0 


2.9 


2.7 





(which is more reliable than plasmon pole models). The 
GW corrections, which differ significantly from what can 
be obtained using a scissor operator, produce BOs in ex- 
cellent agreement with experiment. While the BOs ob- 
tained with and without vertex correction do not differ 
significantly, those resulting from QSGW present a larger 
deviation from experiments. These findings allow us to 
recommend the use of eigenenergy self-consistent GW, 
which is less demanding. Now that the quality of the 
procedure has been assessed, it can be used in a predic- 
tive way for heterojunctions of high technological interest 
such as those in transistors or organic photovoltaics. 
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